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We investigate observational constraints from lensing 
statistics and high-z type la supernovae on flat cosmological 
models with nonrelativistic matter and an exotic fluid with 
equation of state, p x = (m/3 — l)p x - We show that agreement 
with both tests at the 68% confidence level is possible if the 
parameter m is low (m < 0.85) and 0.24 < S7 m o i$ 0.38 with 
lower values of fi m o corresponding to higher m. We find that 
a conventional cosmological constant model with Q m o — 0.33 
is the best fit model of the combined likelihood. 

PACS number(s): 98.80.Hw 



I. INTRODUCTION 

Flat cosmological models with a cosmological con- 
stant are currently the favorite candidates to describe 
the dynamics of the universe. These models are theo- 
retically appealing because A helps to reconcile inflation 
with dynamic estimates for the matter density parame- 
ter (O m o). Further, they fit nicely some observational 
data, as for instance those coming from the classical 
magnitude- redshift test, where Snela are used as stan- 
dard candles, and also from the favored location of the 
first acoustic peak of the cosmic microwave background 
radiation (CMBR) angular power spectra. 

Historically, Einstein himself was the first one to in- 
clude a A-term in the general relativity field equations 
in order to make them compatible with a static universe. 
Several times A was introduced in cosmology but later 
on discarded when improved data became available. The 
preliminary results from Snela almost excluded A, so it 
is curious to note that this time better data is supporting 
A instead of rejecting it From the aesthetic point of 

view, if compared, for instance, to the Einstein-de Sitter 
models, A-models seem ugly. A new parameter is intro- 
duced in the theory and in order to dominate the dynam- 
ics of the universe only in recent times, this parameter 
should have a very small value (A < 10 -56 cm -2 ), which 
is 50 to 120 orders of magnitude below the estimate given 
by quantum field theory. In fact, in part to alleviate this 
problem dynamical-A models were originally suggested. 
There are also observational motivations for considering 
dynamical A models instead of the constant one. For in- 
stance, in these models the COBE-normalized amplitude 
of the mass power spectrum is in general lower than in 
the conventional constant-A model, in accordance with 



observations m. Further, the distance to an object with 
redshift z is smaller than the distance to the same object 
in a constant-A model (assuming the same value of Sl m o). 
So, constraints coming from lensing statistics are weaker 
in these models ppO[|. 

The dynamical-A models present in the literature can 
schematically be divided in three types: scalar field 
x-fluid and decaying- A laws glQ. A 

phcnomenological decaying-A law model in which A de- 
creases as A oc a~ m [here a is the scale factor of the 
Friedman-Robertson- Walker (FRW) metric and to is a 
constant (0 < m < 3)] was suggested in Refs. [pl)|-f2~ij] . It 
was observed that the Einstein equations for these models 
are the same if instead of a A-term, it would be consid- 
ered (beside matter and radiation) a x-fluid with equa- 
tion of state, p x = — l) p x . In spite of the similarity 
at the level of Einstein equations, these two phcnomeno- 
logical models are different. For instance, in the case of 
a decaying A-term, matter is created as a result of the 
decaying vacuum, while in the exotic fluid description 
the x-component is conserved. In this sense the last ap- 
proach is more conventional and closer to the description 
of dynamical-A in terms of a scalar field rolling down a 
potential without coupling to other fields. 

In previous work |pj,|2lj], as a first approximation, 
it was assumed that the x-component is smoothly dis- 
tributed. We were motivated by the fact that if the x- 
component clumps on scales ~ 10 — 20 /i _1 Mpc it would 
be detected by dynamical measurements and this is not 
observed. Another concern is that if the perturbed pres- 
sure is negative (as is the background one), the fluid 
sound velocity would be imaginary, small scales would 
grow exponentially and the system would be highly un- 
stable @|L|,§. However, as observed in Refs. 
the fact that the background pressure is negative not 
necessarily implies imaginary sound speed. Furthermore, 
Caldwell et al. || (see also Ref. [|l0|) pointed out that the 
smoothness assumption is gauge dependent. However, if 
for instance A is modeled by a scalar field, perturbations 
will only give an appreciable effect at large scales. For 
scales well inside the horizon, smoothness is a good ap- 
proximation. 1 16 1 . 

In this paper we shall deal with two cosmological tests: 
gravitational lensing statistics and the Snela magnitude 
redshift test. The strongest observational support for an 
accelerated universe comes from Snela. This test can 
be considered the main motivation for introducing some 
kind of exotic matter with negative pressure. The cos- 
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mological constant is the simplest possibility, but not the 
unique. On the other hand, most lensing statistics anal- 
ysis give lower values for A and we find interesting to 
compare the predictions of these two important tests. 
Here we consider the special case where the exotic com- 
ponent is a x-fluid with constant equation of state and 
that is smooth on scales smaller than horizon. We show 
that agreement with both tests at the 68% confidence 
level is possible if the parameter m is low (m < 0.85) 
and 0.24 < Q m0 < 0.38 with lower values of fl m0 corre- 
sponding to higher m. The best fit model of the com- 
bined likelihood has m — (cosmological constant) and 
fl m0 ~ 0.33 (n A0 ~ 0.67). 

This paper is organized as follows. In Sec. II the basic 
field equations and distance formulae are presented. In 
Sec. Ill we obtain constraints on the models from lensing 
statistics. Constraints from high redshift Snela are ob- 
tained in Sec. IV. We also present a combined likelihood 
analysis of both tests in this section. In Sec. V our main 
conclusions are stressed out. 



II. FIELD EQUATIONS AND DISTANCE 
FORMULAE 

In this paper we consider spatially flat, homogeneous, 
and isotropic cosmologies with nonrelativistic matter and 
an exotic x-fluid with equation of state, p x = (y — l) p x 
(or equivalently a time-dependent A term such that 
A oc a~ m ). Since we are mostly interested on effects 
that occurred at redshift z < 5 we neglect radiation. We 
consider that nonrelativistic matter and x-fluid are sep- 
arately conserved (p m oc a~ 3 and p x oc a~ m ). 

The Einstein equations for the models we are consid- 
ering are: 



fWff 2 (^) 3 + n x0 H 2 ^) m (2.i) 



and 



- _ - tt2 I a 



<^Wfl(£)"\ (2.2) 



where ft m o = 1 — Cl x o is the matter density parameter 
and H is the present value of the Hubble parameter. 

In the next sections we shall use two concepts of cos- 
mological distances, the angular diameter distance and 
the luminosity distance. So, we end this section present- 
ing their definitions and showing some expressions we use 
in our computations. 

Consider that photons are emitted by a source with 
coordinate r = n at time t\ and are received at time to 
by an observer located at coordinate r = 0. The emitted 
radiation will follow null radial geodesies on which 9 and 
<j) are constant. The comoving distance of the source is 
defined by: 



X 



to 



cdt 



(2.3) 



that in flat space is equal to r\ . The present value of the 
scale factor times the comoving distance, gives the proper 
distance, d(0,z), between the source and the observer, 
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and in our case it reduces to the following expression, 

dy 



d(0, z) = cH { 
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The luminosity distance of a light source is defined in 
such a way as to generalize to an expanding and curved 
space the inverse-square law of brightness valid in a static 
Euclidean space, 



d L (0,z) 



C 



1/2 



riao(l + z). 



(2.6) 



Here C is the total energy emitted by the source per 
unit time in its rest frame (absolute luminosity), T is 
the measured flux (energy per unit time per unit area 
measured by a detector). For flat models, d(0, z) = a ri 
and dh can therefore be written as 



d L (0,z) = {l + z)d(0,z). 



(2.7) 



The angular diameter of a light source of proper di- 
ameter D (here supposed to be redshift independent) at 
r = r\ and t = t\, observed at r — and t — to is given 

by 
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The angular diameter distance is defined as the ratio 
of the source diameter to its angular diameter, i.e., it is 
the distance that would be attributed to the light source 
if it were in a Euclidean space, 

d A (0, *) = ? = d(0, z)(l + z)- 1 = d L (0, z)(l + z)- 2 . 
o 

(2.9) 

It is convenient to write the angular diameter distance 
for our models, between the redshift zl and zs, 



d A (z L ,z s ) = 

ell., : '' ; 



1 + z s J ZL 



dy 
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(2.10) 
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III. CONSTRAINTS FROM LENSING 
STATISTICS 

We start defining the following likelihood function, 



Nu 

IB* 



Nr 
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Y[p'ck- 



(3.1) 



k=l 



Here Nl is the number of quasars that have multiple im- 
age, Njj is the number of quasars that don't have, p i <C 1 
is the probability that quasar i is lensed and p ck is the 
configuration probability, that we shall consider as the 
probability that quasar k is lensed with the observed im- 
age separation. To perform the statistical analysis we use 
data from the HST Snapshot survey (498 high luminous 
quasars (HLQ), the Crampton survey (43 HLQ), the Yee 
survey (37 HLQ), the ESO/Liege survey (61 HLQ), The 
HST GO observations (17 HLQ), the CFA survey (102 
HLQ) , and the NOT survey (104 HLQ) We consid- 
ered a total of 862 (z > 1) high luminous optical quasars 
plus 5 lenses. 

The differential probability, dr, that a line of sight in- 
tersects a galaxy at redshift zl in the interval dzL from 
a population with number density uq is, 



dr = 



differential light travel distance cdt 

mean free path 1 / riG~Ka\ 



(3.2) 



where QjQy IS the maximum distance of the lens from the 
optical axes for which multiple images are possible. It 
is a function of the angular diameter distance between 
observer and lens, lens and source, observer and source 
and it also depends on the lens model. 

We use a singular isothermal sphere (SIS) as the lens 
model. Hence, 



cU(0, z L )d A (z L , z s ) 
d A {0,z s ) 



and the bending angle is ao = 4-7r(^) 2 ~ 
1.5 ( 225 ^ m ^ s ) 2 w ith (T|| standing for the one compo- 
nent velocity dispersion. We assume conserved comov- 
ing number density of lenses, rig = no(l + z) 3 and a 
Schechter pfl] form for the galaxy population. So, 



n 



L* y \ L* L* 
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and we take n* = 1.4 ± 0.17 /i 3 10~ 2 Mpc -3 and a = 
— 1.0 ± 0.15 J3^| . The division of the luminosity func- 
tion by galaxy type is taken from Marzke et al. |3lJ , 
n e = 0.61 ± 0.21 /i 3 10~ 2 Mpc~ 3 for early type galaxies 
and n s = 0.79 ± 0.26 /i 3 10- 2 Mpc" 3 for spirals |2f|. As 
a conservative approach biased pro higher values of f^A, 
we do not consider lensing by spiral galaxies. We assume 



that the luminosity satisfies the Faber-Jackson relation 
|§, L/L* = (cr||/(T*)T, with 7 = 4 and take a* = 225 

Km/s. 

The total optical depth (r), obtained by integrating 
dr along the line of sight from to zs, can be expressed 
analytically 

t(z s ) = ^(d A (0,z s )(l + z s )) 3 (cH^)- 3 , (3.4) 

where F = 167r 3 n e (ci/ " 1 ) 3 (cri , j/c) 4 L(l + a + 4/7) ~ 0.026 
measures the effectiveness of the lens in producing mul- 
tiple images |p| . 

It is important to include two corrections to the optical 
depth: magnification bias and selection function due to 
finite resolution and dynamic range p5| . 

Since lensing increase the apparent brightness of a 
quasar and since there are more faint quasars than bright 
ones, there will be over representation of lensed quasars 
in a flux limited sample. The bias factor is given by 

MM 



B(to, z) = Af 2 B(m, z, M , M 2 ) 



(3.5) 



where 



B(m,z,Mi,M 2 ) = 

Since we are modeling the lens by a SIS profile, Mq = 2, 
and we use M 2 — 10 4 in the numerical computation. 

We use the following expression for the quasar lumi- 
nosity function |2q] 



dm V 



(3.7) 



where 



mo + (z + 1) for z < 1, 

777 = { m for 1 < z < 3, (3.8) 

mo — 0.7(z — 3) for z > 3, 

and we assume a = 1.07, b = 0.27 and mo = 18.92. The 
magnification corrected probabilities are 



Pi = T(zi)B(mi,Zi) 



(3.9) 



Finally we have to consider the selection function due 
to finite resolution and dynamic range. It can be shown 
that the selection function corrected probabilities are: 

n 

p , K z) = ft J # g( mi ^p a ) (3 10) 



and 



Pci =Pci{0) 



B(m,z,M ,M 2 ) 

pA B(m,z,M f (8),M 2 ) 
p'J B(m,z,M ,M 2 ) 



(3.11) 
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where ff54[ 
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M/(0) = M 
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/> 1 



A, (3.12) 



(3.13) 



and 



/ = f(6) = io°- 4Am W. 



(3.14) 



To simplify computation we use two selection functions, 
one for the HST observations and another one for all the 



ground based surveys 33 1 . Using more accurate selection 
functions for each ground based observations separately 
have little statistical effect. 

Recently Falco et al. J2^] observed that statistical lens- 
ing analysis based on optical and radio observations can 
be reconciled if the existence of dust in E/SO galaxies is 
considered. In our computation we assume a mean ex- 
tinction of Am =0.5 mag as suggested by their estimates. 

By expressing Ci ens as a function of the parameters 
to and £l m o we obtained the maximum of the likelihood 
function (£Hs) and formed the ratio I = C iena It 
can be shown that with two parameters, the distribution 
of — 21n^ tends to a x 2 distribution with two degrees 
of freedom p8| , p5] |. In two dimensions the 68%(lcr) and 
95.4% (2a) confidence levels are where the likelihood is 
60.7% and 13.5%, respectively, of the peak likelihood. 

In Fig.l we plot contours of constant likelihood (95.4% 
and 68%) in the two parameter space (m, f2 m o). The 
maximum of the likelihood occurs for m ~ 2.4 and £l m o = 
0. The same approach when applied to constant A models 
(since to = we now have only one degree of freedom) 
gives: fl A < 0.76 (or fi ro0 > 0.24) at 2cr, 1 > Q m o £ 0.39 
at lcr with a best fit at O m o ~ 0.62. For comparison, 
without considering extinction, in the case to = 0, we 
get: Qa ^ 0.55 at 2a, that is slightly more conservative 
than Kochanek's ]25| fl\ < 0.66 at 2a. 

In Fig. 2 contours of constant number of multiple im- 
ages as a function of the parameters to and Q m o are dis- 
played for the case in which Am = 0.5 mag extinction is 
considered. Comparing the two figures we see that the 
la contour of Fig. 1 corresponds roughly to the contour 
in Fig. 2 where 9 multiple images are expected. For the 



best fit model in two dimensions the predicted number 
of multiple images is ~ 5.3. 

The results obtained in this section are more accurate 
than those presented in Ref. 0] . We now took into ac- 
count magnification bias in the configuration probabil- 
ity and considered the selection function due to finite 
resolution and dynamic range. In the present approach 
extinction is also considered, so the constraints are less 
restrictive with respect to A than those obtained in Ref. 
M ■ 
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FIG. 1. Contours of constant likelihood (95.4% and 68%) 
arising from lensing statistics. 
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FIG. 2. Contours of constant number of multiple images as 
a function of the parameter space. 



IV. CONSTRAINTS FROM HIGH-REDSHIFT 
TYPE IA SUPERNOVAE 

There are two major ongoing programs to systemati- 
cally search and study high-z supernovae. Although the 
very preliminary results indicated a low value for the 
cosmological constant ( 57a < 0.51 at the 95% confi- 
dence level) jj), more recent analysis with larger sam- 
ple of supernovae, now points to a different direction. 
Now the data indicate an accelerated expansion such that 
51a ~ 0.7, f2 m o ~ 0.3 and strongly supports a flat Uni- 
verse. Hi]. 

In this section we update the constraints on the equa- 
tion of state obtained by Silveira & Waga (|lj , where the 
first data from the Supernovae Cosmology Project was 
used S. We now consider data from the High-z Super- 
novae Search Team. We use the 27 low-z and 10 high-z 
Snela (we include SN97ck) reported in Riess et al. jjj and 
consider data with the MLCS J|6|,[l| method applied to 
the supernovae light curves. The results for spatially flat 
models that we present here, are similar to those obtained 
by Garnavich et at. j35fl . 

Following a procedure similar to that described in 
Riess et al. |lj] , we determine the cosmological parame- 
ters m and fl m .o through a x 2 minimization neglecting 
the unphysical region f2 m o < 0. To simplify computation 
we fix the Hubble parameter to Hq = 65.2 km/s Mpc -1 
[ 1 1 , but the results are independent of this choice for Hq 
[l,p5|. We use 
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where 



[i, p = 5\ogd L + 25, 



(4.1) 



(4.2) 



is the distance modulus predicted by each model, /j,o is 
the observed (after corrections) distance modulus, a^ Q its 
uncertainty and a v is the dispersion in galaxy redshift 
due to peculiar velocities. Following |l]] we use a vi = 
_5 200km/s jjjp-jj.g g ne j a w ith z not derived from 

In 10 czi & 

emission lines in the host galaxy we add 2500 km/s in 
quadrature to 200 km/s (see Table 1 in [0). 
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FIG. 3. Contours of constant likelihood (95.4% and 68%) 
arising from lensing statistics (dashed lines) and type la su- 
pernovae are shown. 

In Fig. 3 contours of constant likelihood 95.4% (2a) and 
68% (la) arising from the Axlne analysis are displayed 
together with those from lensing (dashed lines). For 
Snela the peak of the likelihood is located at to ~ 1.1 and 
Q m0 = 0. If we fix to = we get O m0 = 0.25 ± 0.08 (la) 
(for comparison Riess et al. obtained Q m o — 0.24 ±0.1). 
From the figure it is clear that there is a region in the 
parameter space (the region inside the triangle with ver- 
tices (to ~ 0.85,fi m0 ~ 0.24), (to = 0,n ro0 - 0.32) and 
(to = 0,n m o — 0.38)) such that all points are inside the 
la (68%) confidence region of both tests. 

In Fig. 4 we display contours (95.4% and 68%) of the 
combined (lensing plus Snela) likelihood. For the com- 
bined x 2 analysis we used Xtot = ^Xsne — 21nZ, with 
I = Ciens/C^ns as defined in Sec. III. Although the 
peak of the likelihood for each test separately occurs at 
!l m o = 0, the maximum of the combined likelihood oc- 
curs at m — (cosmological constant) and £l m o — 0.33. 
Note that best fit models of the combined likelihood are 
in accelerated expansion (go < 0). Models with m = 2 
(cosmic strings ]12[|) and any value of £l m o are at more 
than 99% c.l. away from the peak of the likelihood. 
We observe that if, for instance, we take h = 0.65 and 
n B h 2 = 0.02, the CMBR first acoustic peak (£ poa k), for 
models with to and fl m o inside the la allowed region in 
Fig 4, will have £ U cak values between ~ 215 and 230, (see 
Fig. 4 in Ref. [fig] ), that are close to the current best 
values for £ pca ^ obtained from CMBR data Models 
with parameters to and f2 m o in this region are in agree- 
ment with the current CMBR data as well. Constraints 
from observations of clusters suggest fi m o > 0.14 J37|. In 
Fig. 4 we display this constrain as a dotted line. Models 
above this line are preferred. 
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FIG. 4. Contours of combined likelihood (95.4% and 68%) 
arising from lensing statistics and type la supernovae are 
shown. 



V. SUMMARY 

We have studied observational constraints from lens- 
ing statistics and high-z Snela on spatially flat cos- 
mological models whose matter content is nonrelativis- 
tic matter plus an exotic fluid with equation of state, 
Px — (^p — Using a lensing approach, where ex- 

tinction is considered, and that takes into account mag- 
nification bias and the selection function due to finite 
resolution and dynamic range in the configuration prob- 
ability, we obtained the 68% and 95.4% confidence level 
on the parameters m and fi TO o. The present results up- 
date those obtained by Silveira & Waga pi]] . We also 
update the Snela constraints on the equation of state 
obtained in Ref. where the first data from the Su- 
pernovae Cosmology Project was used. We considered 
data from the High-z Supernovae Search Team, the 27 
low-z and 10 high-z Snela reported in Ref. jjj. We used 
data with the MLCS method applied to the supernovae 
light curves. We showed that agreement with both tests 
at the 68% confidence level is possible if the parameter 
m is low (m < 0.85) and 0.24 < £l m o ^ 0.38 with lower 
values of Q m o corresponding to higher m. We observed 
that best fit models of the combined likelihood are in 
accelerated expansion (go < 0). We obtained that a con- 
ventional cosmological constant model with Q m o a 0.33 
(f^Ao = 0.67) is the best fit model of the combined likeli- 
hood. 
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